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Spin- flip hot spots, i.e., special k-points on the Fermi surface showing a high spin-mixing param- 
eter, are found on the surface Brillouin zone boundary of ultrathin films of monovalent (noble and 
alkali) metals, in contrast to the long-known fact that they do not occur in the bulk of these metals. 
Density-functional calculations within the Korringa-Kohn-Rostoker Green function method in Cu, 
Ag, and Au 10-layer (OOl)-films show that the hot spots can have a substantial contribution, e.g. 
49% in Au(OOl), to the integrated spin-mixing parameter, that could lead to a significant enhance- 
ment of the spin-relaxation rate or spin-Hall angle. Owing to the appearance of spin-flip hot spots, a 
gigantic anisotropy of the Elliott- Yafet parameter [52% for Au(OOl)] is also found in these systems. 
The findings are important for spintronics applications in which the dimensionality of spin probes 
is reduced. 

PACS numbers: 72.25.Rb, 73.50.Bk, 72.25.Ba, 85.75.-d 



The term of spin-flip hot spots was coined by Fabian 
and Das SarmsP^ when, based on the Elliott-Yafet 
mechanismj^ they predicted a high spin- flip probability 
for electrons at certain special points on the Fermi sur- 
face of bulk Aluminum. The hot spots are formed in the 
vicinity of band degeneracies or near-degeneracies, fre- 
quently occuring at the Brillouin zone boundaries, high- 
symmetry lines or accidental degeneracy points. On the 
other hand, hot spots should be absent in monovalent 
metals,!^ because their Fermi surface is nearly spherical 
and in most cases does not cross the Brillouin zone edge. 

These conclusions were derived for the bulk.^ ^ ^ How- 
ever, in the case of ultrathin films spin-flip hot spots 
may occur even in monovalent metals owing to the dif- 
ferent shape of the surface Brillouin zone compared to 
the bulk. This is the main conclusion of the present pa- 
per, which we base on simple theoretical arguments and 
demonstrate by ah initio calculations in Cu, Ag, and Au 
ultrathin films in the (001) crystalline orientation. We 
also investigate the contribution of the hot spots to the 
Elliott-Yafet parametei^l^ I? (defined below) as well as 
to its anisotropy. Our findings should be accounted for in 
spintronics applications where the hot spots play a role, 
such as giant magnetoresistance, spin Hall effect and spin 
dynamics or even quantum computation.l^^KI^ 

We proceed with a short summary of the theoreti- 
cal background.^ ^ ^ In non-magnetic systems with space- 
inversion symmetry and in the presence of the spin-orbit 
coupling the Bloch states at any k-point are at least 
twofold degenerate and comprise a superposition of spin- 
up and spin-down states that is frequently called spin 
mixing: 

*+(r)= [ak(r)|t)+5k(r)|;)]e*-- 

(1) 

*k(r)= [alk(r) l;) - 61k(r) |t>] e^'^'-. 

Here, the spinors |t) and \\) are eigenvectors of the 
z-component Pauli matrix (given a chosen z-axis) 
and ak(r) and 6k (r) are the lattice-periodic parts of the 



Bloch states corresponding to the large and small spin- 
component, respectively. The spin-expectation value 
of these partner states is S(k) = (^klf^l^k) = 
~ (^k I f ^l^k ) with a the vector of Pauli matrices. The 
index refers to the state with maximal z-component 
of the spin-expectation value ^'^(k) = ~ ^k)' where 
we defined the space-integrated spin-mixing parameter 
^k = / <^^^|6k(r)P G [0,^]. Maximizing the value of 
^'^(k) is done with respect to all possible linear combina- 
tions of the two degenerate states at k, and it is obvious 
that one can choose a different pair of states in Eq. ([T]) 
by maximizing the projection of S(k) along any spin- 
quantiation axis s prescribed by the experimental con- 
ditions (e.g. by an external magnetic field or b y the po- 
larization direction of an injected spin current ).'m^ We 
will return to this freedom of choice below when defining 
the anisotropy. The relation between the large and small 
components of and ^^j^ implied in Eq. ([I]) follows 
from time-reversal and space-inversion symmetry.^ ^ It is 
also convenient to define the Elliott-Yafet parameter as 
the Fermi-surface average: 

where dS\^ is the Fermi-surface element, is the Fermi 
velocity, n{EY) is the density of states at the Fermi level 
Ey and V^z is the Brillouin zone volume. As the value 
of iP' depends on the choice of s along which S(k) is 
maximized, we define the anisotropy of the Elliott-Yafet 
parameter as^^ 

^ _ maxs 6^(s) - ming 6^(s) 
mius 6^(s) 

The spin-mixing parameter is the main quantity 
of interest in the analysis of many spin-flip related phe- 
nomena, as it reflects the deviation of a Bloch state from 
being a spin eigenstate. When 6^ happens to be large 
(close to ^), then the particular state has an almost 
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completely mixed spin character Concerning spin re- 
laxation, an electron scattered into this state practically 
loses its spin character — for example, according to El- 
liott's approximation^ the spin-relaxation time Ti is in- 
versely proportional to the Elliott-Yafet parameter 6^. 
Concerning the spin-Hall effect, such a state with large 6^ 
exhibits a strong contribution to the conductivity tensor 
due to the high value of the Berry curvature associated 
with itPi^ The k-points where b^. becomes large are the 
spin-flip hot spots that we are examining in the present 
work. The fact that these occur at degeneracy or near- 
degeneracy points, in particular at the back- folding of 
the energy bands at the Brillouin zone edge, follows from 
a consideration of the transitions from a band at energy 
to a band at + Ak for small inter-band separation 
Ak under the action of the spin-orbit Hamiltonianl^J^ 
From these arguments it also follows that monovalent 
metals should not show spin- flip hot spots, since their 
almost spherical Fermi surface either does not cross the 
Brillouin zone boundary (as in all alkali metals except 
Cs) or, if it crosses the boundary (as in the noble metals, 
forming a neck around the L-point), then it does so in a 
way that there occur no hot spots. ^ 

However, the situation of monovalent metals changes 
when one considers ultrathin films. In this case the two- 
dimensional periodicity implies a surface Brillouin zone, 
while the Fermi surface consists of rings occuring when 
the spherical Fermi surface of the bulk system is cut par- 
allel to the surface plane at positions determined by the 
film-induced quantization in the perpendicular direction. 
This well-known effect is schematically demonstrated in 
Fig.[l^. Considering a (001) film, where the surface Bril- 
louin zone is a square, the larger of Fermi rings exit 
the first Brillouin zone and a back-folding occurs (see 
the dashed lines and their back-folded counterparts in 
Fig. This, under action of the periodic potential, 

can form a weak degeneracy lifting at the Brillouin zone 
edge, with the resulting states being energetically very 
close. But this is precisely the case when spin-flip hot 
spots occur under the action of the spin-orbit coupling 
between the near-degenerate states. 

To what extent this actually happens in a realistic 
case depends on the exact shape of the Fermi surface, 
the number of film layers, the surface orientation, and 
of course the material. Here we present ab initio re- 
sults on Cu, Ag and Au (001) films of 10 layers thick- 
ness where we find an effect of considerable magnitude, 
while in alkali (001) thin films the magnitude of the ef- 
fect depends very sensitively on the film thickness due 
to the weak spin-orbit coupling. The electronic struc- 
ture is calculated within the local density approximation 
to density-functional theory.l^ For our calculations we 
employ the full-potential Korringa-Kohn-Rostoker Green 
function method using the experimental lattice parame- 
ter and ignoring surface relaxation; we use an angular 
momentum cutoff of /max = 3. Details on the formalism 
and implementation can be found in Refs. [T5HT7I 

Starting our analysis with the bulk properties, the 




FIG. 1: Sketch of the ideal Fermi surface of a monovalent 
metal, (a) Spherical in bulk; the circles shown at certain 
values of kz correspond to the quantization in a film grown in 
the ^-direction, (b) Forming circles in a (001) film and back- 
folding in the surface Brillouin zone, shown together with the 
crossed surface Brillouin zone boundary. 



spin-mixing parameter for bulk Cu and Au is well stud- 
ied e.g. in Ref. 18 and Ref. |T5| It is long-known that 
the Fermi surfaces of the noble metals are similar. The 
texture of 6^ on the Fermi sufaces is also similar, but 
the magnitude is very different due to the much stronger 
spin-orbit coupling in Au. The value of the spin-mixing 
parameter on the Fermi surface in bulk Cu varies between 
0.0005 and 0.002, in Ag between 0.0008 and 0.0025, while 
in Au it clearly reaches much higher values varying be- 
tween 0.01 and 0.045. It is also found that there are no 
hot spots on the noble-metal Fermi surfaces. The Fermi- 
surface average is calculated to be 6^(Cu) = 0.0015, 
b^{Ag) = 0.0017, and ^^(Au) = 0.03. 

Now we examine the distribution of b^ on film Fermi 
surfaces. We should point out from the beginning that 
the values of 6^, as well as the integrated 6^, depend on 
the direction of the spin-quantization axis s, because the 
matrix elements of the spin-flip part of the spin-orbit op- 
erator between Bloch states change with respect to the 
s. This anisotropy effect, its microscopic origin and its 
relation to hot spots was analyzed in Ref. 11 for bulk ma- 
terials, but it is certainly present also for films [see Ref.fTOl 
for an analysis in W(OOl)]. Practically in the case of films 
the effect is that, if the spin-quantization axis s is chosen 
in-plane ([100]-direction), the value of b^ is different com- 
pared to the case when the procedure is carried out with 
the spin-quantization axis chosen perpendicular to the 
film along [001]. Physical consequences of the anisotropy 
are, e.g., a variation of the spin-relaxation time or of 
the spin-Hall conductivity with respect to the direction 
of polarization of the spin current in the material, since 
the polarization direction corresponds to s in the present 
theory. In bulk Cu, Ag, and Au, the cubic symmetry 
and the absence of hot spots makes the anisotropy of 6^ 
negligible (less than 0.1% comparing s ||[001], [110] and 
[111]) but in the films it takes gigantic values due to the 
hot spots. 

In Fig. [2] (top and middle), the distribution of spin- 
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spin-flip hot spots 
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FIG. 2: (Color online) Top: Distribution of spin-mixing pa- 
rameter on the Fermi surface of a Au(OOl) 10-layer film with 
the spin-quantization axis along [001], i.e. out of plane (left 
panel) and in-plane along the [100] direction (right panel). 
The full surface Brillouin zone is shown. Middle: Focus on 
the Brillouin zone edge of Au(OOl) shown in the top panels 
in order to distinguish the extent of the hot spots. An asym- 
metry of distribution of bl- can be seen in the case of s || [100] 
(right). Bottom: Distribution of spin-mixing parameter on 
the Fermi surface of a Cu(OOl) (left) and a Ag(OOl) (right) 
10-layer film with the s || [001] (out of plane). 



mixing parameter on the Fermi surface of a 10-layer 
Au(OOl) film with the spin-quantization axis along [001] 
(left) and along the [100]-direction (right) is shown in a 
color code. It is easily seen (see middle panel for a fo- 
cused view) that spin-flip hot spots are found near the 
edge of Brillouin zone. While for most k-points inside 
the Brillouin zone the spin-mixing parameter has a value 
less than 0.05, at the edge of the Brillouin zone it ex- 
hibits very high values reaching even the maximal value 
of fully spin-mixed states, b^. = ^. The anisotropy is 
already exposed by the observation that the color-coded 
texture of bl in the left panel of Fig. [2] (s || [001]) is dif- 
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TABLE I: Values of the Elliott- Yafet parameter (xlO^) in 
bulk and in 10-layer Cu, Ag and Au (OOl)-films. In the case 
of the films the values are given for three directions of s with 
respect to the crystallographic axes ([001] corresponds to the 
direction of s perpendicular to the film). The anisotropy val- 
ues A correspond to the films. The anisotropy in bulk is less 
than 0.1%. 



ferent (and more symmetric) than the one on the right 
panel (s || [100]), especially if one focuses on the region 
close to the hot spots. 

Taking the average over the Fermi surface, we obtain 
62([001]) = 0.064 for Au(OOl) with s || [001]. This value 
is considerably higher than the value of 0.03 that we 
find for bulk Au or of 0.035 that we find for Au(lll) 
10-layer films; the latter show no hot-spots but have 
Rashba-type surface states^^ with an enhanced spin- 
mixing parameter. On the contrary, in the (001) thin 
films, there are no surface states at the Fermi level, yet 
the value of b^ is higher. This goes to show the impor- 
tance of the spin-flip hot spots in this case. To estimate 
the contribution of spin-flip hot spots to the total spin- 
mixing parameter, we perform the integration in Eq. ([2| 
only for the k-points belonging to a small region around 
the hot spots. The results show that 49% of spin-mixing 
parameter in Au(OOl) comes from the hot spots. 

In the case that s is along the [100] direction (Fig. [2j 
top and middle panels, right) the values of b^ change; 
the maximum value is then 0.44 and the integrated one 
is 6^ ([100]) = 0.042. We obtain for the anisotropy A = 
[62([001]) -62([100])] /62([100]) a value of 52%. This 
value is gigantic compared to the negligible anisotropy 
in the bulk of noble metals, and is comparable in mag- 
nitude to the anisotropy in e.g. W(OOl) films^^ where it 
arises from surface states or to the anisotropy in the bulk 
of hep Os (59%)^^ where it arises from larger spin- flip hot 
areas. 

For Cu and Ag (001) thin films, the spin- flip hot spots 
are also present at a similar position near the surface Bril- 
louin zone edge as in Au (the Fermi surfaces of Cu(OOl) 
and Ag(OOl) 10-layer films are shown in Fig. |2j bottom 
panel). The hot spots are, however, less intense due to 
the weaker spin-orbit coupling of Cu and Ag. The con- 
tribution of these spin-flip hot spots to the total 6^ is 
nearly 10%, which is again smaller than Au(OOl) case. 
The anisotropy A of the Elliott-Yafet parameter is found 
to be 7% and 37% of Ag(OOl) and Cu(OOl) 10-layer films, 
respectively. The Fermi surface averaged 6^ in Cu(OOl) 
is 0.002 which is also somewhat larger than the value of 
0.0016 found in Cu(lll) 10-layer thin films. Our results 
are summarized in Table [H 
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Concluding the discussion on the noble metal films, we 
comment on the absence of hot spots in the (111) film 
orientation. Part of the reason is that the fcc(lll) atomic 
planes are more close-packed than the (001) atomic 
planes, resulting in a somewhat larger surface Brillouin 
zone. Another part of the reason is that the (111) sur- 
face Brillouin zone is hexagonal, and in this sence closer 
in shape to the maximal circle forming the equator of the 
bulk Fermi surface. As a result of both, the projection 
of the bulk Fermi surface almost fits in the surface Bril- 
louin zone, leaving only little room for crossing the zone 
boundary. If such crossings appear and lead to hot spots 
has to be tested in each material and thickness separately, 
but as we find, it is not the case here. 

Finally, for completeness, we briefly discuss the spin- 
mixing enhancement in alkali- metal thin films, although 
they are typically not used in spintronics devices. Even 
for the alkali metals the Fermi surface has a non- 
vanishing p and d character that is responsible for spin- 
orbit coupling with strength = {l\ 2m?c^r ^\^) (^^^^ 
|/) the p {I = 1) or d {I = 2) wavefunction and V the 
potential). However, the nuclear potential (that causes 
most of the spin-orbit coupling) is well screened by the 
filled shells of the core electrons, contrary to the noble 
metals where a larger part of the screening is done by the 
valence electrons and by the not-fully-localized d band. 
Additionally, the d character of the noble-metal Fermi 
surfaces is more pronounced compared to the alkali met- 
als and the spin-orbit coupling of the d states is stronger 
in the noble metals because of the high localization of the 
d bands. As a consequence, the spin-orbit strength in al- 
kali metals is significantly lower than in the noble metals 
and the hot spots are much less pronounced. Still, we 
found that at some film thicknesses, e.g. 10 layers of 
Rb(OOl), it so happens that the Fermi surface without 
spin-orbit coupling is degenerate at the Brillouin zone 
edge meaning that the first Fourier component of the pe- 



riodic potential vanishes (at least to numerical accuracy, 
where we cross-checked this result using the full-potential 
linearized augmented plane wave method^^). In this case 
spin-orbit coupling causes a splitting with full spin mix- 
ing when s is perpendicular to the film, i.e. hot spots 
with hl^ = \ emerge at the Brillouin zone edge. Yet the 
magnitude of 6^ drops very quickly as the bands separate 
with increasing distance from the edge, e.g. 6^ = 0.02 al- 
ready at a distance of 0.0005 x (where aiat is the 
lattice constant). As a result, the integrated parameter 
iP' remains very low. 

In summary, we have shown that the Fermi surface 
of monovalent metals in an ultrathin film geometry can 
show spin-flip hot spots as the Fermi rings cross the sur- 
face Brillouin zone boundary. This i s in c ontrast to the 
bulk of such metals, where it is knowri^^^ that hot spots 
do not occur, as the Fermi surface is included within 
the Brillouin zone. We have furthermore shown that 
the spin-flip hot spots contribute to gigantic anisotropy 
values of the spin-mixing parameter with respect to the 
relative orientation between the spin-quantization axis 
and the crystallographic directions. Since the presence 
of hot spots strongly influences the spin-relaxation time 
or the spin-Hall conductivity, our findings can have con- 
sequences in spintronics applications, in particular if ul- 
trathin noble-metal films are used to transmit or probe 
spin currents. The calculated anisotropy can very likely 
lead to a variation of the spin-relaxation time or spin- 
Hall conductivity with respect to the spin-polarization 
direction of the spin current in experiments. 
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